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Abstract 

A search for the decay K® — > \i~ is performed, based on a data sample of l.Ofb -1 
of pp collisions at \fs = 7 TeV collected by the LHCb experiment at the Large Hadron 
Collider. The observed number of candidates is consistent with the background-only 
hypothesis, yielding an upper limit of B(Kg — > < 11(9) x 1CT 9 at 95 (90)% 

confidence level. This limit is a factor of thirty below the previous measurement. 
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1 Introduction 



The decay K[ 



is a Flavour Changing Neutral Current (FCNC) transition that 



has not yet been observed. This decay is suppressed in the Standard Model (SM), with an 
expected branching fraction [Tl|2] 



b(k: 



s ° ^ /iV) = (5.0 ±1.5) x 10" 
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while the current experimental upper limit is 3.2 x 10 7 at 90% confidence level (C.L. 

Although the dimuon decay of the K® meson is known to 
B(K° n+fT) = (6.84 ±0.11) x 10" 9 @, in 
new particles can still be observed in K® — > \i 
decay width of K® s — > can be written as [5] 



El- 
be 

agreement with the SM, effects of 
/i" decays. In the most general case, the 



o 

L,S 



8tt 



2m^ 
m K 



\A\ 2 +\l 




(1) 



where A is an S-wave amplitude and B a P-wave amplitude. These two amplitudes 
have opposite CP eigenvalues, and in absence of CP violation (K® = K®, K® = K®), 
K® decays would be generated only by A while K® decays would be generated only by 
B. The decay width T(K® — > /i + /U~) receives long-distance^ contributions to A from 
intermediate two-photon states, as well as short distance contributions to the real part 
of A. In any model with the same basis of effective FCNC operators as the SM, the 
contributions from B can be neglected for B(K® — > fi + fi~). The decay width of K® — > /x 4 fi~ 
depends on the imaginary part of the short-distance contributions to A and on the long- 
distance contributions to B generated by intermediate two-photon states. Therefore, the 

fi~ ) in agreement with the SM does not necessarily imply that 
has to agree with the SM. Contributions up to one order of magnitude 
above the SM expectation are allowed [2j ; enhancements of the branching fraction above 
10~ 10 are less likely. The study of K® — > has been suggested as a possible way to 

look for new light scalars p]. 

In addition, bounds on the upper limit of B(K® — > fi + fi~) close to 10~ n could be very 
useful to discriminate among scenarios beyond the SM if other modes, such as K + — > -k + vv, 
were to indicate a non-standard enhancement of the s — > dll transition^] 

The KLOE collaboration has searched for the related decay — > e + e~, which is 
affected by a larger helicity suppression than the muonic mode, and set an upper limit on 
the branching fraction B(K® — > e + e~ 



measurement of B(K® — > fi 



B(K° S /i 



< 9 x 10~ 9 at 90% confidence level |6|. 
The LHC produces ~ 10 13 per fb _1 inside the LHCb acceptance. In this paper, 
a search for K® — > ji + /i~ is presented using 1.0 fb _1 of pp collisions at y/s = 7 TeV 
collected by LHCb in 2011. Dimuon candidates are classified in bins of a multivariate 



1 The long-distance scales correspond to masses below that of the c quark, while short-distance scales 
correspond to masses of the c quark and above. 

2 Charge conjugation is implied throughout this paper. 
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discriminant, and compared to background and signal expectations. The background 
present in the signal region is a combination of combinatorial background and K® — > tt + tt~ 
decays in which both pions are misidentified as muons. The number of expected signal 
candidates for a given branching fraction hypothesis is obtained by normalising to the 
measured — > 7r + 7r~ rate. The results obtained by the measurements in different bins 
are combined, and a limit is set using the CL S method (7^8). The data in the signal region 
were only analysed once the full analysis strategy was defined, including the selection, the 
binning and the evaluation of systematic uncertainties. 

The LHCb apparatus, and the aspects of the trigger relevant for this analysis are 
presented in Sect. |2j Section [3] is devoted to the full signal selection and to the definition of 
the multivariate method used as the main discriminant. In Sect. |4]the different backgrounds 
for K® — > decay are described, as well as the expected background in the signal 

region. The normalisation, required to convert the number of K® — > n + fi~ candidates to 
the branching fraction, is detailed in Sect. |5j The systematic uncertainties are described 
in Sect. [6} The limit setting procedure, together with the corresponding expected and 
observed limits, is presented in Sect. [7], and conclusions are drawn in Sect. [8j 

2 Experimental setup 

The LHCb detector [9] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
(VELO) surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes placed downstream. The combined tracking 
system has a momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 
100 GeV/ c, and an impact parameter (IP) resolution of 20 |xm for tracks with high transverse 
momentum (pp) with respect to the beam direction. Charged hadrons are identified using 
two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter and 
muon systems, followed by a software stage which applies a full event reconstruction. For 
this analysis, the events are first required to pass a hardware trigger which selects muons 
with > 1.5GeV/c. In the subsequent software trigger [10] , at least one of the final state 
tracks is required to be of good quality and to have > 1.3 GeV/c, an IP > 0.5 mm and 
the x 2 °f the impact parameter (IP x 2 ), defined as the difference between the \ 2 °f the 
primary vertex built with and without the considered track, IP x 2 > 200, where the impact 
parameter is calculated with respect to the associated primary vertex in the event. Events 
with high multiplicity are rejected using a series of global event cuts. A prescale factor 
of 0.5 is applied to the lines triggered by the K® — > fi + fi~ candidates. The K® — > 
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candidates responsible for the trigger of both the hardware and software levels are called 
TOS (trigger on signal). Notice that the LHCb trigger has not been optimised for the 
K® — > search, but rather to the decays of the heavier c and b hadrons. 

Events with a reconstructed K® — > fi + fi~ candidate can also be triggered independently 
of the signal candidate if some other combination of particles in the underlying event 
passes the trigger. Such candidates are called TIS (trigger independently of signal). 
The TIS and TOS categories are not exclusive as muons from both the Kg — > 
candidates and from the underlying event can pass the trigger. There is overlap between 



the two, which allows the determination of trigger efficiencies from the data 11 . Finally 



a Kg — > /i yU~ candidate may be triggered by a random trigger; such candidates are called 
minimum bias (MB). As the sample produced by this trigger is negligible in terms of a 
potential Kg — > signal, this is only used to provide Kg — > 7r + 7r~ candidates, useful 

to understand the distributions that the signal would have in the case of no trigger bias. 

The K^ — > 7r + 7r~ decay is used as a signal proxy in this analysis, with the exception 
of the muon particle identification. The Kg — > 7t + tt~ and Kg — > fi + fi~ candidates are 
reconstructed requiring two opposite charge tracks with hits in both the VELO and 
tracking stations. Approximately 40% of the Kg mesons inside the LHCb acceptance decay 
within the VELO. Additionally, for Kg — > / u + pr candidates, hits in the muon chambers 
are required. The LHCb mass resolution for the combination of such tracks (ranging from 
3.2 to 5MeV/c 2 depending on their kinematics and very similar both for Kg — > tt + tt~ 
and Kg — > fi + n~) allows signal to be discriminated from Kg — > 7t + tt~ decays where both 
pions are misidentified as muons. The Kg — > 7t + tt~ (in the hypothesis) and the 

Kg — > mass peaks are separated by around 40MeV/c 2 (see Fig. [I]). 

The signal candidates are selected requiring a decay vertex (secondary vertex, SV) 
well separated from the pp interaction point (primary vertex, PV) with a Kg candidate 
reconstructed momentum pointing to the PV, and an invariant mass in the region of the 
Kg mass. 

For the simulation, pp collisions are generated using Pythia 6.4 [12] with a specific 
LHCb configuration [13] . Decays of hadronic particles are described by EvtGen [14] 



in which final state radiation is generated using Photos 15 . The interaction of the 



generated particles with the detector and its response are implemented using the Geant4 



toolkit 16 as described in Ref. 171. 



3 Selection and multivariate classifier 

A set of requirements is applied on the reconstructed candidates to reduce the size of the 
sample to be analysed. This selection is identical for the signal and for the Kg — > 7r + 7r~ 
control channel, with the exception of the particle identification requirements for the 
daughter tracks and the mass window for the Kg. The selection includes requirements on 
the quality of the daughter tracks (x 2 < 5 per degree of freedom), on their IP \ 2 > 100; 
and on the distance of closest approach between them (< 0.3 mm). Requirements on the 
decay time (> 8.9 ps, with the SV being downstream of the PV) and on the IP of the Kg 
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Figure 1: Mass spectrum for selected K® — > 7r + 7r~ candidates in the MB sample. The 
black points correspond to the mass reconstructed under the tttt mass hypothesis for the 
daughters, while the red triangles correspond to the mass reconstructed under the \i\x 
mass hypothesis. 



meson (< 400 um) are also applied. Furthermore, A decays are vetoed via a requirement 
in the Armenteros-Podolanski plane |18|, by including cuts in the transverse momentum of 
the daughters with respect to the K® flight direction and in the longitudinal momentum 
asymmetry of the daughters. Figure [l] shows the mass spectrum for selected K® — > 7r + 7r~ 
candidates in the MB sample after applying the set of cuts described above. 

In order to further increase the background rejection, a boosted decision tree (BDT) [19 
with the AdaBoost algorithm |20| is used. To avoid possible biases (in particular, as 
K® candidates in the mass sidebands are also used to get the background predictions) 
the sample is split in two subsamples, for which the BDT discriminants are optimised 
separately. The BDT trained in the first half-sample is used as a discriminant for the 
second half, and vice versa. This guarantees the absence of an overlap between the samples 
used for the training and those used in the analysis, while making maximal use of the 
available data both for the BDT training and for the background evaluation. The BDT 
discriminants are trained on data, using K® mass sidebands as a background training 
sample and K® — > 7t + tt~ decay as a signal training sample. In addition, these samples are 
further split according to the TIS and TOS trigger categories in order to take account of the 
bias introduced by the trigger on the variable distributions. Thus, in total, four different 
samples are defined (two subsamples for the TIS trigger category and two subsamples for 
the TOS category) and are combined as described in Sect. [7j The candidates that are 
simultaneously TIS and TOS are analysed only as TIS candidates to avoid counting them 
twice. Only one per mille of the TOS candidates overlap with TIS candidates. In order to 
minimise the difference in the momentum distributions of the daughters of the candidates 
from the TOS, TIS, and MB samples, and hence to ease the normalisation to K® — > 7r + 7r~ 
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decays, for the TOS analysis at least one of the daughters is required to have a px above 
1.3 GeV/c. 

The variables entering in the BDT discriminants are: 

• the decay time of the K® candidate, computed using the distance between the SV 
and the PVj^] and the reconstructed momentum of the K® candidate; 

• the smallest muon IP \ 2 of the two daughter tracks with respect to any of the PVs 
reconstructed in the event; 

• the Kg IP x 2 with respect to the PV; 

• the distance of closest approach between the two daughter tracks; 

• the secondary vertex x 2 , which adds complementary information with respect to the 
distance of closest approach of the tracks, as it uses information on the uncertainty 
of the vertex fit; 

• the angle of the decay plane in the K® rest frame with respect to the K® flight 
direction, which is isotropic for signal decays, but not necessarily for background 
candidates; 

• variables used to discriminate against material interactions, as further detailed below. 

An important source of background consists of muons resulting from interactions 
between the particles produced in the PV and the detector material in the region of the 
VELO. The position of the SV of the background candidates from the K® mass sidebands 
in the x — z plane is shown in Fig. [2] The structures observed correspond to the position 
of the material inside the VELO detector. To discriminate against this background, two 
different approaches are used for the TIS and TOS categories, consisting of two different 
choices of variables for the BDT. 

For the TOS category, two additional variables are included in the BDT, the px of the 
K® and a boolean matter veto that uses the VELO geometry to assess whether a given 
decay vertex coincides with a point in the detector material or not. Muons from material 
interactions have a harder px spectrum than muons from other background sources and 
hence are more likely to be selected by the trigger. The use of this variable in the BDT 
gives slightly better performance than simply applying the veto as a selection. 

For the TIS category, the coordinates of the position of the SV in the laboratory rest 
frame are used to deal with the background from material interactions. As the simultaneous 
use of the lifetime, px of the K% meson, and the SV position allows the BDT to effectively 
compute the mass of the candidate, a fake signal peak could be artificially created out 
of the combinatorial background. Hence the px of the K® meson is not used in the TIS 
analysis. This second approach is found to be better than the matter veto (and K® px) 
for the TIS analysis, while, on the contrary, the matter veto boolean variable is found to 
be better than the SV coordinates for the TOS analysis. 

3 If more than one PV is reconstructed, the PV associated to a given particle is the one that minimises 
the IP x 2 of that particle. 
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Figure 2: Position in the x — z plane of the secondary vertices of the background candidates 
found in the high mass sideband for (left) TIS candidates and (right) TOS candidates. 
The lighter coloured areas correspond to higher density of points. 

The BDT discriminants for the TIS category are trained on the K® — > 7r + 7r~ TIS 
sample, which is used as a proxy for signal decays, and on the upper mass sideband of TIS 
candidates as proxy for the background candidates. The region above 1.1 GeV/c 2 (above 
the resonance) is used to decide on the BDT settings and the region between 504 MeV/ c 2 
and 1.0 GeV/c 2 (below the resonance) is used to train the BDT algorithm chosen. The 
K® — > 7r + 7r _ decays with pions misidentified as muons and passing the same trigger as 
Kg — > signal are used as a proxy for the signal decays for TOS candidates. In order 

to minimise differences between misidentified K® — > 7r + 7r~ decays and true K® — > 
decays, tight muon identification requirements are applied to the K® — > 7t + tt~ misidentified 
sample. These tight requirements are chosen such that the trigger biases in the K® — > tt + tt~ 
simulated decays were the same as in the K® — > simulated decays. The upper mass 

sideband of TOS candidates is used as a proxy for background candidates with the same 
splitting of the mass regions used for the TIS category. 

Candidates with low values of BDT response are not considered because of the large 
amount of background at low BDT values. This BDT cut causes about 50% efficiency loss, 
depending on the sample. The rest of the candidates are classified in ten bins of equal 
signal efficiency, i.e. a total of forty bins are combined to get the CL a limit. 

4 Background 

The search region is defined as the mass range [492, 504] MeV/c 2 . The background level is 
calibrated by interpolating the observed yield from mass sidebands ([470, 492] and [504, 600] 
MeV/c 2 ) to the signal region. This is done by means of an unbinned maximum likelihood 
fit in the sidebands, using a model with two components. The first component is a power 
law that describes the tail of K® — > 7t + tt~ decays where both pions are misidentified as 
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Figure 3: Background model fitted to the data separated along (left) TIS and (right) TOS 
trigger categories. The vertical lines delimit the search window. 



muons; this model has been checked to be appropriate using MC simulation. The second 
component is an exponential function describing the combinatorial background. As an 
illustration, Fig. [3] shows the distribution of candidates for all BDT bins and for TIS and 
TOS samples, respectively. 

Other sources of background, such as K® — > 7r + /i~P M , K® — > /i + /i~7, K® — > /i + /i~7, 
K® — > 7r + /i - z/ M and K® — > fi + fi~ decays, are negligible for the current analysis. The 
contamination from K® — > decays may be relevant for future analyses as it can 

produce an effective K® — > branching fraction up to ~ 2 x 10 -11 . 



5 Normalisation 

A normalisation is required to translate the number of K® — > / u + signal decays into a 
branching fraction measurement. Two normalisations are determined independently for 
TIS and TOS candidates. The B(K® — > is computed using 

B{K° S -> fi + fj-) = e wn N K o^ +tl - 

where, in a given BDT bin, N K o_+ + _ is the observed number of signal decays, iV is: o_ >7r + 7r - 
the number of — > 7r + 7r~ decays, and e TO /e w the ratio of the corresponding efficiencies. 
The efficiencies are factorised as e = e SEL e PID e TRIG / SEL where: 

• e SEL is the offline selection efficiency. It includes the geometrical acceptance, recon- 
struction and selection, i.e, it is the probability for a K® — > 7r + 7r~ (K® — > ji + //~) 
decay generated in a pp collision, to have been reconstructed and selected; 

• e PID is the efficiency of the muon identification for reconstructed and selected 
K® — > fi + fi~ signal decays; 
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e TRic/SEL = ^sel&pid&trig/^sel&pid where t/RIG denotes either the TIS or the 
TOS categories, is the trigger efficiency for decays that would be offline selected. 
Under this definition, trigger efficiencies can be determined from data using the 



procedure described in Ref. 11 



The ratio of reconstruction and selection efficiencies between K® — y fi + n~ decays and 
Kg — y 7r + 7r~ decays is evaluated using simulation. In order to account for data/simulation 
differences in the differential K® cross-section, this efficiency is evaluated in bins of pt and 
rapidity of the K®. The reconstruction and selection efficiency for K® — y tt + tt~ decays 
is between 60% and 85% (depending on which point in the phase space a given event is 
from) of that of the K® — y fi + fi~ decays due to difference in the material interactions of 
the pions compared to muons. 

The factor e PID is evaluated in bins of the BDT (both for the TOS and TIS categories) 
by measuring the muon identification efficiency as a function of p and px using calibration 
muons. The sample of calibration muons is obtained from a J ftp — y sample in which 

positive muon identification is required for only one of the tracks. The p and p^ spectra of 
the pions from K® — y tt + tt~ decays in a MB sample is later used to get the efficiency for 

— y fi + fi~ decays. The e PID efficiency is between 68% and 82% (depending on the BDT 
bin and the sample). It is measured with a precision between 1% and 10%. For the ratio 
of trigger efficiencies, different strategies are considered for the TIS and TOS samples. 

For the TIS samples, the K® —y ji + \i~ yield is normalised to the K® — y TIS 
yield. In this case, the trigger efficiencies cancel in the ratio, because the probability to 
trigger on the underlying event is independent of the decay mode of the K® meson. This 
cancellation is verified in simulation. The normalisation expression for TIS decays reads 

n(K Q ^n+u-\ ^SEL i TIS/SELjyTIS 

B(K» -> 71+71-) ~ C SEL e PID TIS/SEL jyTIS > W 

where g TIS / SEL denotes the trigger efficiency for signal and K® — y 7r + 7r~ modes respectively 
and N T1S is the number of TIS decays in a given BDT bin for signal and K® — y 7r + 7r~ 
modes respectively. 

For the TOS sample, the K® —y fi + n~ yield is normalised to the K® —y 7r + 7r~ yield 
from MB triggers. To perform such a normalisation, the global event cuts (the trigger 
requirements which reject high multiplicity events) that affect the lines triggered by the 
Kg — y candidates are applied offline to the MB events, so that their effect cancels. 

The normalisation requires in this case an absolute determination of the TOS trigger 
efficiency for K® — » fi + /j,~ , e™ S/,SEL , as well as the knowledge of the average prescale factor 
of the MB trigger, s MB . The absolute TOS trigger efficiency for the signal, in bins of BDT, 
is computed using muons from B + — y J/ip(-^y ^i + ^)K + decays]^] The p and p^ spectra of 
the B + — y J/ip{— y ^i + ^r)K + muons are reweighted in order to match those of pions from 
the K® — y 7r + 7r~ decays. Trigger unbiased p and spectra of the K® — y 7r + 7r~ decays can 



4 To avoid bias, it is required that another object be the origin of the trigger and not the muons alone, 
i.e. the muons from this sample are TIS. 
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be obtained from the MB sample. The normalisation expression for TOS decays reads 



g(gg-^+;Q _e^ 1 , MB 

B(K° -+ TT+TT") £ SEL £ PID TOS/SEL tf**^ > ^ > 

N^_^ n+7T _ being the number of — > tt + tt~ decays from the MB trigger and N™!^, + - 
denoting the number of signal decays from the TOS category. 
The quantities 

TIS f SEL ,PID TIS/SEL /V TIS ^ ' 

MM MM ^MM ifg->7r+7r- 



and 



e SEL ! .MB g(g^^-) 
"TOS e g EL e p ID TOS/SEL jyMB \°) 



are called normalisation factors and are defined for each of the BDT bins. For a given 
number N of —> /i + ' \i~ signal decays, the corresponding value of B[K\ —> is 
then a x N. Using the value of B(K® — > 7r + 7r~) from Ref. |4j, the normalisation factors 
are in the range [6.6, 16.2] x 1(T 8 for the TIS category, and [0.9, 7.8] x 1(T 8 for the TOS 
category, depending on the BDT bin. 



6 Systematic uncertainties 

The quantities considered in the determination of the branching fraction that are affected 
by systematic uncertainties are listed below. 

• The background expectations per bin, obtained by comparing the results with the 
model described in Sect. [3] to those computed: a) if the combinatorial background is 
modelled by a linear function; b) if the mass range over which the fit is performed 
is modified; c) repeating the fit excluding (together with the signal region) the 
12MeV/c 2 left and right windows neighbouring the search window and comparing 
the fit prediction to the yields in those regions; no correlation is considered among 
the different bins for this systematic uncertainty. 

• The ratios of reconstruction and selection efficiencies and absolute muon identification 
efficiencies, for which systematic uncertainties are obtained from the difference 
between different methods and from the comparison to simulation (around 20% for 
the ratios and 5% for muon identification efficiencies); no correlation is considered 
among the different bins. 

• The branching fraction of the normalisation channel B(K® — > 7r + 7r~) = (69.20 ± 
0.05)% [4 1; its uncertainty affects coherently the signal expectations of the forty bins 
of the analysis. 
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• The absolute TOS efficiency, for which the systematic uncertainty is obtained from 
the comparison to simulation (around 15%); no correlation is considered among the 
different bins. 

• The effective prescale factor of the MB sample, s MB = (2.70 ± 0.76) x 10" 6 . The 
uncertainty is evaluated from the difference between the prescale factor as measured 
in data and the value of the prescale as set in the trigger system. This systematic 
uncertainty affects coherently the signal expectations of the twenty bins of the TOS 
analysis. 

The leading systematic uncertainties are those coming from the ratio of trigger efficien- 
cies for the TOS analysis and from the ratio of reconstruction and selection efficiencies for 
the TIS analysis. 



7 Results 

The modified frequentist approach (or CL S method) (Tj,[8] is used to assess the compatibility 
of the observation with expectations as a function of B(K® — > fi + fi~). 

Test-statistics are built from pseudo-experiments for the signal plus background and 
background-only hypotheses. For each pseudo-experiment a product of likelihood ratios is 
computed depending on the expected number of signal events for a given branching fraction, 
Si, the expected number of background events, bi and the observed number of events, 
di for bin i. The CL s+b (CL b ) is defined as the probability for signal plus background 
(background only) generated pseudo-experiments to have a test-statistic value larger than 
or equal to that observed in the data. The CL S is defined as the ratio of confidence levels 



CL 



s+b 



CL 



b 



This ratio is used to set the exclusion (upper) limit on the branching fraction, 
whereas 1 — CL b is used as a p- value to claim evidence or observation. A 95(90)% confidence 
level exclusion corresponds to CL S = 0.05(0.1). 

The values of bi are obtained from the fit of the mass sidebands, as detailed in Sect. |3J 
The values of Sj depend on the assumed branching fraction, as well as on the normalisation 
factors computed in Sect. [5j The uncertainties on the input parameters are taken into 
account by fluctuating the signal and background expectations when generating the b 
and s + b ensembles. These fluctuations are performed via asymmetric Gaussian priors, 
following the formula 

x'i = Xi ^1 + ^r(s + - s_) + \r 2 (s + + s_)^ (7) 
where Xi is the central value of the parameter, r is a random number generated from a 



normal distribution and s + and s_ are the relative (signed) errors of Xi 21 . Correlations 
are implemented by using the same value of r for the parameters that should fluctuate 
coherently. 

The observed distribution of events is compatible with background expectations, giving 
a p-value of 27%. Figure [|] shows the expected and observed CL S curves for the TIS 
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Figure 4: CL S curves for (a) TIS, (b) TOS categories and for (c) the combined sample. 
The solid line corresponds to the observed CL S . The dashed line corresponds to the median 
of the CL S for an ensemble of background-alone experiments. In each plot, two bands 
are shown. The green (dark) band covers 68% (la) of the CL S curves obtained in the 
background only pseudo-experiments, while the yellow (light) band covers 95% (2<r). 



category and for the TOS category as well as for the combined measurement. The upper 
limit found is 11 (9)xl0 -9 at 95 (90)% confidence level and is a factor of thirty below the 
previous world best limit. Table [T] summarises the limits in the TIS, TOS categories, and 
the combined result. 



8 Conclusions 

A search for K® — > has been performed in l.Ofb -1 of data collected at the LHCb 

experiment in 2011. The candidates observed are consistent with the expected background, 
and the upper limit 

B(K° /i+/0 < 11(9) x 10~ 9 
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Table 1: Upper limits on B(K® — > n~) for the TIS and the TOS categories separately, 
and for the combined analysis. The last entry in the table is the p- value of the background- 
only hypothesis. 



Quantity TIS TOS Combined 

Expected upper limit at 95 (90)% C.L. [10~ 9 ] 42 (33) 13 (10) 11 (9) 

Observed upper limit at 95 (90)% C.L. [10~ 9 ] 24 (19) 15 (12) 11 (9) 

p-value 0.95 0.20 0.27 



at 95(90)% confidence level is an improvement of a factor of thirty below the previous 
world best limit [3]. The observed p- value for the background only hypothesis is 27%. 
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